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O ■ ABSTRACT 

liHj' We study the Faraday rotation measure (RM) due to the Galactic magnetic 

I ■ field (GMF) toward high Galactic latitudes. The RM arises from the global, reg- 

■ ular component as well as from the turbulent, random component of the GMF. 
Q>^ ■ We model the former based on observations and the latter using the data of 
*^ ■ magnetohydrodynamic turbulence simulations. For a large number of different 
^ ■ GMF models, we produce mock RM maps around the Galactic poles and calcu- 
O ■ late various statistical quantities with the RM maps. We find that the observed 
^ . medians of RMs toward the north and south Galactic poles, ~ 0.0 ±0.5 rad m~^ 

■ and ~ +6.3 ± 0.5 rad m~^, are difficult to explain with any of our many alter- 
. nate GMF models. The standard deviation of observed RMs, ~ 9 rad m~^, is 

^ . clearly larger than that of simulated RMs. The second-order structure function 

■ ■ ■ of observed RMs is substantially larger than that of simulated RMs, especially 

at small angular scales. We discuss other possible contributions to RM toward 
high Galactic latitudes. Besides observational errors and the intrinsic RM of 
background radio sources against which RM is observed, we suggest that the RM 
due to the intergalactic magnetic field may account for a substantial fraction of 
the observed RM. Finally we note that reproducing the observed medians may 
require additional components or/and structures of the GMF that are not present 
in our models. 
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Introduction 



Cosmic magnetism, the origin and nature of magnetic fields in ou r universe, is one of 
outstanding problems of modern astrophysics (see iGaensler et al.ll2004[ ). Exploration of the 
Galactic magnetic field (GMF) and the intergalactic magnetic field (IGMF) is listed as one 
of key science projects for the Square Kilometer Array (SKA), and is one of the important 
science projects for Jansky Very Large Array (JVLA), Murchison Widefield Array (MWA), 
the Low Frequency Array (LOFAR), the Australian SKA Pathfinder (ASKAP), and the 
South African Karoo Array Telescope (MeerK AT) (see, e.g.. ICarilli fc Rawlingsll2004j : iBeck 



20091 : iKrause et al.l l2009l : iGaensler et al.l 120101 ). Measuring Faraday rotation, the rotation 



of the plane of linearly-polarized radio emission due to the birefringence of magneto-ionic 
medium, is a powerful method to observe the GMF as well as the IGMF. Up to now, there 
have been a number of studies of Faraday rotation measure (RM) to elucidate the structure 
and statistical properties of the GMF and the IGMF. 



For the IGMF, there h ave been a number of observations of RM in clusters of galaxies 



see 



Carilli &: Taylorl l2002l for a review). Observations have been also extended to the 



outside of clusters, toward the cosmic web, and extr agalactic contr i butions of RM due to 



the IGMF or others have been discussed (see, e.g.. IXu et al.l 120061 : iKronberg et al.l 12008 



filarn e nts of galaxies (Rvu et al. 



2010l : lAkahori fc Rvu 



2010 



Il998a 



Bernet et al.l l2012l : iHammond et al.l 120121). There are seyeral theoretical works for RM in 



Dolag et al. 



201 ll ). For instance. 



2005 



Cho fc Rvu I2OO9I: IStasyszyn et al. 



Akahori fc Rvul J201oh inv estigated RM 



in filaments using a model IGMF based on a turbulent dynamo (iRyu et al.l l2008l ). They 
found that the root-mean-square (rms) valu e of RMs through a sin gle filament is expected 
to be ~ 1 rad m~^ in the local universe. lAkahori Sz Ryu! (|201l[ ) extended this work by 
using the redshift distribution of polarized background radio sources against which RM is 
observed. They found that the rms value of RMs through filaments up to redshift ~ 5 would 
be ~ several rad m~^. They also found that the second-order structure function (SF) has 
a nearly-fiat profile in angular separations of > 0.2°, meaning that RMs through filaments 
decorrelate on angles less than ~ 0.2°. 

Recently, the RMs toward the North and South Galactic poles (N GP and SGP, r espec- 
tively) have been investigated to study the GMF as well as the IGMF. iTaylor et al.l (|2009[ ) 
studied the RM data from the NRAO VLA Sky Survey (NVSS), and estimated non-zero 
vertical strengths of the GMF, about -0.14 ±0.02 /iG and +0.3 ±0.03 /xG toward the NGP 
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and SGP, respectively. iMao et al.l (120 lOf ) used RM data from the Westerbork Radio Syn- 
thesis Telescope (WSRT) and the Australia Telescope Compact Array (ATCA), and found 
that the median value of RMs toward the SGP is +6.3 ± 0.5 rad m^^ (corresponding to a 
vertical GMF strength of +0.31 ±0.02 fiG), while that toward the NGP is 0.0 + 0.5 rad m^^ 
(+0.00 + 0.02 /iG). The standard devi ations of RMs we re ~ 9.2 rad m~^ and 8.8 rad m~^ to- 



ward the NGP and SGP, respectively. iMao et al.l (120101 ) also put an upper li r nit of ~ 1 fiG on 
the strength of random magnetic field at high Galactic latitudes. IStil et al.l (120111 ) examined 
the NVSS data in detail, and found that the second-order SFs of RM at angular separations 
of > 1° have a value ~ 100-200 rad^ m"^ toward the NGP and ~ 300-400 rad^ m"^ toward 



the SGP. In addition, based on the latitude dependence of RM, ISchnitzelerl (120101 ) examined 
the Galactic and extragalactic contributions to RM in the NVSS data. He estimated that 
the Galactic contribution is o"rm,mw ~ 6.8 + 0.1(8.4 + 0.1) rad m~^ and the extragalactic 
contribution (including those intrinsic to the polarized background radio sources and due to 
the IGMF) is ctrm,eg ~ 6.5 + 0.1(5.9 + 0.2) rad m~^ for the northern (southern) hemisphere. 

The separation of the Galactic and extragalactic contributions in observed RMs, how- 
ever, is not trivial. It requires good understandings of the GMF and its contribution to ob- 
served RMs. There have been a number of theoretical works to model the Ga l actic RM (e.g.. 



Sun et al. 



2011 



20081: 



Wael 



Pshirkov et al 



cens et all bood: ISun fc ReichI bood: Ijaffe et al.l I2OI0I: IVan Eck et aL 



20 



cosrn i c-ray propagati on: 



iMao et al 



2012 



Jansson fc Farrai 



201^ 



(see also works on 



Prouza fc Smidal l2003i r iTinyakov fc TkachevI l2005l : iGiacinti et al.l 



2OIOI : iTakami &: Satdl2010l ). The Galactic RM arises from the global, regular component as 
well as from the turbulent, random component of the GMF. While the regular component 
has been modeled with analytic fitting formulae based on observations, the random compo- 
nent h as been modeled usin g power-law spectra with random phases in Fouri er space. For in 



stance, 



Sun fc ReichI (120091 ) used the pubhcly-available HAMMURABI code (IWaelkens et al. 



20091 ). and adopted a Kolmogorov-like power spectrum with average amplitude 3 fiG in a 
box of 10 pc size. They found that at Galactic latitudes |6| ~ 70°, the second-order SF has 
a magnitude of up to a few xlOO rad^ m~^ at angular scales of > 10'. 

Previous studies have successfully reproduced the observed properties of the Galactic 
RM as well as those of the radio continuum emission toward low and mid-Galactic latitudes 



(e.g.. lSun et al.ll2008l ). There is, however, a lack of studies that can be compared with recent 
observations of RMs toward the Galactic poles. In this paper, we simulate the Galactic RM 
toward high Galactic latitudes and investigate its statistical properties. While we model the 
regular component of the GMF based on a number of observations, we use the data of three- 
dimensional magnetohydrodynamic (MHD) turbulence simulations to model the random 
component. We obtain simulated maps of the Galactic RM toward the Galactic poles in 
a field-of-view (FOV) of 900 deg^ to compare with recent observations and in a 200 deg^ 
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FOV to predict future SKA observations. We estimate how much the GMF can contribute 
to observed RMs and discuss how the statistics of simulated RMs compare with those of 
observed RMs. The rest of the paper is organized as follows. In Section 2, we describe our 
models. The results are shown in Section 3. Discussion and summary follow in Sections 4 
and 5, respectively. 



Models 



2.1. Regular Components 

We first briefly describe our models for the global, regular components of the electron 
density and the GMF. The details along with the definitions of the coordinate systems used 
in the paper are described in Appendices. We also briefly describe models for the electron 
temperature and the rms speed of random flow motions, which are used in modeling the 
electron density fluctuations and the turbulent, random magnetic field in the next subsection. 



For the electron density, n^fi, we employ the NE2001 model (ICordes fc Lazidl2002U2003[ ) 
displayed in Figure [TJ We take the original parameters of the model except for the scale 
height of the thick disk, hi, and the mi d-plane electron dens ity for the thick disk, rieX, we 
use hi = 1.8 kpc and rig,! = 0.014 cm~^ (iGaensler et al.l 120081 ) . which better reproduce both 
the dispersion measure (DM) and emission measure (EM) toward high Galactic latitude in 
our model (Appendices). The top- left panel of Figure [2] shows one-dimensional profiles of 
ne,o from the Sun along the NGP and SGP. The electron density in the northern sky is 
smaller than that in the southern sky at low altitudes. This is due to the presence of a 
l ocal, hot bubble with ne,o = 0.005 cm~^, centered at {x,y,z) = (0.01,8.45,0.17) in kpc 
(ICordes fc L"^l2002h . 



For the regular magnetic field, Bq, we use combination s of the a x i-syrn metric spi- 
ral (ASS) or bi -symmetric s p iral ( BSS) field introduced by ISun et al.l (120081) . th e halo 
to roidal field of ISuii &: ReichI (|2010[ ) , and the halo poloidal fields of iGiacinti et al.l (120 lOf ) 



or 



Jansson &: Farrarl (120121 ) . The top panels of Figure [3] show one-dimensional profiles of the 
strength of the regular field, Bq, and the line-of-sight (LOS) field strength, from the Sun 
along the Galactic poles. For instance, ADPN indicates a model toward the NGP (N) in- 
cluding the ASS field (A) and the dipole toroidal field (D) with the dipole poloidal field (P). 
For the full list of models and the model name convention, refer to Table 1. We note that 
the dipole and quadrupole models are identical toward the NGP, but differ toward the SGP. 
The DS models tend to have larger Bq than other models, since in other models the spiral 
and toroidal fields partly cancel each other. The spiral field dominates Bq at \z\ < 1.25 kpc. 
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while the toroidal field dominates at \z\ > 1.25 kpc. For the LOS field strength, the field 
strength is B\\ ~ 0.2 — 0.3 fiG if the poloidal field exists. Otherwise, B\\ = 0. 

For the electro n temperature distribution, we employ the analytic expression adopted 
bv lSun et all J2008h 

Te(i?, z) = 5780 + 287i? - 526|z| + 17702^ (1) 

where Tg is in units of K and the Galactocentric cylindrical coordinates R and z are in kpc. 
The bottom-left panel of Figure [2] shows the one-dimensional profile from the Sun along the 
Galactic poles. 

For the rms speed of random fiow motions, Kms? Ha o bservations provide a weak con- 
straint on the plausible range of values, Kms ^ 15 — 5 km/s (ITufte et al 



1999; Haffner et al. 



20031 : iHill et al.ll2008l : iHaffner et al.ll2010h . iHill et all ( 120081 ) studied the distribution of emis- 
sion measure (EM), and found that the rms Mach number Mj-ms = K-ms/c^ ~ 1.4 — 2.4 for 
|6| > 10° and Mrms is smaller at higher Galactic latitudes. Studies of polarization gradients 
also broadly constrain Mrms ~ 0.5 — 2 (IGaensler et al.ll2011t iBurkhart et al.ll2012l ). We em- 
ploy the simplest uniform distribution of V^ms = 15 km/s (Mrms ~ 0.2 — 1) or 
Vrms = 30 km/s (Mfins ~ 0.5 — 2), as shown in the right panels of Figure [2J The bottom 
panels of Figure |3] show one-dimensional profiles of (3o, the ratio of the gas pressure to the 
magnetic pressure due to Bq. Our f3o is in the range of ~ 0.01 — 100. Here, we assume that 
hydrogen is fully- ionized, helium is neutral, and that their mass fractions are X = 0.76 and 
Y = 0.24, respectively. 



2.2. Random Components 



Observations suggest the presence of electron density fiuctuations and turbulent mag- 
netic fields, in addition to smooth components in the Galaxy. For instance, the volume 
filling factor of electrons, ~ 0.05 — 0.5, which quantifies the dum piness, has been estimated 



from dispersion rneasur e s (DMs) and einissions/absorptions ( e .g.. IPeterson fc Webberl 12002 



Berkhuiisen et al.l 120061 : iHill et al.l 120081 : IGaensler et al.l 120081 : ISun et al.l 120081 ). Highly dis- 



turbed dist ributions of RM and polarization angle clearly indicate turbulent structures of t he 
GMF (e.g., ISun et al.ll2008l : iTavlor et al.ll2009[ IWaelkens et al.ll2009t lOaensler et aPboilh . 



The random, turbulent components of the electron density and the GMF could be 
modeled analytically w ith preassigned spe ctra and random phases in Fourier space. As 
mentioned in Sectio n l.lSun &: Reich (|2009|) followed such an appr oach with a constant am- 
plitude everywhere. iJaffe et al.l ( 120101 ) and lJansson fc Farrarl ( 120121 ) introduced an "ordered" 
or "striated" field component. These works successfully reproduced the Galactic RM mainly 
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toward low and mid-Galactic latitudes. However, it is not clear whether such treatments 
are good enough for studies of the Galactic RM toward high Galactic latitudes, where the 
random components are the dominant contribution to the RM. And constant rms ampli- 
tudes for electron density fluctuations and turbulent magnetic field would not be justified 
for broad distributions of Mrms and /3o toward high altitudes (see Figures |2] and [3]) , since 
they depend on Mrms and /3o (see Appendix). In addition, in turbulent flows, phases are not 
really random. 

An alternative approach would be to use the random electron density given by the 
NE2001 model for the electron density fluctuations. But we would then need to model the 
turbulent magnetic field separately. We therefore do not adopt this approach. 



As in lHill et al.l (|2008[ ). we model the random components with MHD turbulence simula- 
tions in a closed box, as described below. In principle, if we performed full MHD simulations 
of the Galaxy, we could reproduce the electron density fluctuations and the turbulent mag- 
netic field as well as the regular components in the Galactic disk and halo. But currently 
available computational resources do not allow a numerical resolution high enough to simul- 
taneously reproduce both the large-scale global components and the small-scale turbulent 
components. 



2.2.1. MHD Turbulence Simulations 



We embed the data of MHD turbulence simulations in the Galactic halo described by the 
global distributions of riefl and Bq of the previous subsection. For this purpose, we carried out 
three-dimensional simulations of driven, isothermal, compressible MHD turbulence without 
self-gravity, usin g a multi-dimens ional MHD code based on the Total Variation Diminishing 
(TVD) scheme (IKim et al.lll999l ). This is an MH D extens i on of the explicit, second-order 
finite-difference, upwinded, conservative scheme of iHartenl (119831 ) for hydrodynamics. The 
version of the code used i ncludes a flux con straining scheme that maintains V ■ -B = up 
to the machine accuracy (iRyu et al.lll998bl ). In uniform, static medium with Bq assumed 
to be along the x-direction, turbulence was driven by imposing a solenoidal forcing with 
perturbations drawn from a Gaussian random fleld at wavenumbers around fcdrive in Fourier 
space. The amplitude of the forcing was flxed in such a way that turbulence saturates at 
preassigned Mach numbers. Simulations were done in a periodic computational box with 
512^ grid zones. 

The size of the simulation box, Lbox, and the outer scale or the driving scale of tur- 
bulence, Ldrive = 27r//cdrive, Were determined as follows. Observations suggest that Ldrive is 
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on the order of ~ 1 p c and ~ 100 pc in spiral arms and in interarm regions, respectively 



(iHaverkorn et al.ll2008l ). These are ascribed to stellar sources such as stellar winds and pro- 
tostellar outflows, or supernova and superbubble explosions. ISun &: ReichI (120091 ) adopted 
Lhox = 10 pc for studies of RM in low and mid- Galact ic latitudes (I/drive is not defined in 
their approach). On the other hand, iHill et al.l (120081 ) studied the EM distribution esti- 
mated with the Wisconsin Ha Mapper (WHAM). Using MHD turbulence simulations, they 
found that the EM distribution depends on Lbox and L^me] larger Lbox and Ldrive results 
in smaller EMs. They argued that the histogram of the WHAM data on the warm ionized 
medium that includes high Galactic latitude data could be reprodu ced with Lbox ~ 500 pc 
and Ldrive ~ 50 — 250 pc. Adopting the results of iHill et al.l (120081 ) . we set Lbox = 500 pc 
and Ldrive — -^box/2 = 250 pc. The grid of 512^ uniform zones for the box of Lbox = 500 pc 
corresponds to the spatial resolution of ~ 1 pc. We note that this resolution is enough to 
resolve the characteristic scales of turbulence; for instance, the most ener gy containing scal e 
and the integral scale, LkE{k) ~ 50 pc and Lint ~ 75 pc, respectively (see lCho fc Ryull2009l ). 
are sufficiently large compared to the spatial resolution. 

Representative simulations were performed for M^^as = 0.5, 1, 2 and /3o = 0.1, 1, 3, 10, in 
order to cover the ranges of Mrms and /3o toward the NGP and SGP (Figures [2] and [3]). Table 
2 lists the simulations along with the rms values of density and magnetic field strength. For 
the case of (3o < 0.1 at high Galactic altitudes {\z\ > 5 kpc. Figure [3]), we use the data for 
Po = 0.1. Since the contribution from high altitudes to the integrated RM is at most several 
per cent, this approximation does not seriously affect our results. In turbulence si mulations, 
V„ns saturates roughly in one flow-crossing time, ~ Lbox/K-ms (IWu et al.l |2012| ). For the 



length scale of 500 pc, the corresponding timescale of t ~ 1.6 x 10^ yr and ~ 0.8 x 10^ yr 
for the rms flow speed of 15 km/s and 30 km/s is sufficiently short compared to the age of 
the Galaxy. We thus use the data at the saturation stage of turbulence. 



2.2.2. Construction of Model Space by Data Stacking 

The data of simulations are stacked (or piled up) from the Galactic mid-plane up to the 
outer edge of the NE2001 model, \z\ = 10.0 kpc (Figured]). In the stacking, we try to align 
the direction of the regular field in simulation data with that of the analytic model as follows. 
As we already noted, the regular field is dominated by the spiral field at \z\ < 1.25 kpc and 
by the toroidal field at \z\ > 1.25 kpc (see Figure [3]). In each domain, the regular field 
direction is close to either the +x or —x direction in our coordinate systems, except narrow 
transition regions with a width ~ 0.2 kpc. 
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We divide each of the northern and southern hemispheres into four blocks: 



block 1 


(0.00 < 


\z\ 


< 1.25), 


block 2 


(1.25 < 


\z\ 


< 2.50), 


block 3 


(2.50 < 


\z\ 


< 5.00), 


block 4 


(5.00 < 


\z\ 


< 10.0), 



where z is in units of kpc. We calculate the average values of Mj-ms and (5q in each block and 
seek the simulation of closest parameters. Table 1 summarizes the direction of the regular 
magnetic field and the values of Mrms and /3o adopted in each block. For example, in block 
1 of AD0N15, where 15 means Kms = 15 km/s, we choose the simulation with Mrms = 1-0 
and /So = 3.0, and align the regular field with the +a;-direction. 

The stacked data are renormalized in such a way that we match the average electron 
density and magnetic field in the computational box with those of the regular components 
of Section 2.1. After the renormalization, each grid keeps the fluctuations based on MHD 
turbulence simulations. 



2.2.3. Goodness of the Modeling of Random Components 



Examples of profiles of the resulting electron density, magnetic field, and RM (see Section 
2.3) from the Sun along the NGP, including both the regular and random components, are 
shown in Figure IH It is apparent that fluctuations in the case with K-ms = 30 km/s are larger 
than those for V„ns = 15 km/s, as expected. The strength of the turbulent magnetic field 
{h = \B — Bq\) is at most a few /xG, and mostly < 1 /xG. This seems to be smaller than the 
strerigth of the randorn field, a few to several /xG, in the Galactic disk (e.g., lOhno &: Shibata 
19931 : iBeck et al.l Il996[ ). but consistent with ~ 2 /xG a few kpc from the Galactic plane 
estimated by assuming equipartition between the thermal and nonthermal pressures ( ICox 
20051 ). as well as w ith the recent es timates of < 1.5 /iG and < 1.4 /iG toward the NGP and 
SGP, respectively (IMao et al.ll2010( ). The sign of LOS magnetic field changes several times. 
|RM| increases with the path length, and then saturates around \z\ ~ 2.0 kpc where the 
electron density becomes small. 

We check our model for the electron density distribut i on by comparing the resulting 



DM a nd EM with observed values (IPeterson k Webberll2002l : iHill et al.ll2008l : iGaensler et al. 



20081 ). DM, which is determined only by the global, regular component, is 22 and 25 pc cm 



toward the NGP and SGP, respectively, while the observed value is ~ 23 — 26 pc cm ^ . 
We note that the NE2001 model, from which we get the regular component, was designed 



- 9 - 



to reproduce observations including DM. The EM in our model is 1.5 (1.0) and 2.3 (1.6) 
pc cm~^ toward the NGP and SGP, respectively, for V].^s = 30 (15) km/s, which is in a 
good agreement with the observed value, ~ 1 — 2 pc cm~^. Ou r mod el also reproduces the 
overall distribution of EM at \b\ > 60° presented by iHill et al.l (120081 ). Details are given in 
the Appendix. 



2.3. Rotation Measure Map 



The RM toward a source outside the Galaxy is defined as 

RM (rad m'^) = 0.81 / n^B\\ds, 



(3) 



where Ue, B 



and the path length s are in units of cm~^, /iG, and pc, respectively. Con- 
ventionally RM is positive when the LOS magnetic field points toward us, and Smax is the 
maximum distance along the LOS up to the outer edge of our model space. In this paper, 
we calculate the RM only due to the GMF. We do not include other possible contributions, 
such as the RM due to the IGMF, the intrinsic RM at sources, and observational errors. We 
discuss those in Section 4. 

We consider a square FOV of 900 deg^ with 30° on a side centered on each of the 
Galac tic poles, which i s comparable to the FOVs of RM observations toward the NGP and 
SGP (IMao et al.l l2010l : IStil et al.l 1201 ll ). In addition, we also consider a square FOV of 
200 deg^ with 14.14° on a s ide. The smaller FO V is roughly that proposed for the dense 
aperture array of the SKA (IFaulkner et al.ll2010l ). The smaller FOV is also used to check 
the FOV dependence of RM statistics. To see structures i n angular separations of > 0.3° fo r 
which the second-order SF of observed RMs is available (IMao et al.ll2010l : IStil et al.l 1201 ll ). 
we choose the angular resolution of pixels to be < 0.3°. The numbers of pixels we use are 



pix 



256^ and 128^ for 900 deg^ and 200 deg^ FOVs, respectively. 



We calculate RMs with Equation 02]) assuming one extragalactic radio source in each 
pixel, and construct mock RM maps. For each map, we randomly rotate stacked simulation 
boxes around the axis parallel to the regular magnetic field in Table 1 (keeping the direction 
of the regular magnetic field) and randomly shift box centers to avoid the repeat of the same 
grid zones which would make artifacts in statistics. When the integration reaches the top of 
a simulation box, the integration proceeds into the next stacked box and the integration is 
continued. If the integration reaches the side of the simulation box, we apply the periodic 
boundary condition to the side, and replicate the box beyond the side. 



We obtain 200 mock maps for each model listed in Table 1 and each FOV, and calculate 



statistical quantities. The quantities shown below in Sections 3.3 and 3.4 are the averages 
for 200 RM maps. 



3. Results 

3.1. Two-Dimensional Map 

Figure O shows sample maps over a 30° x 30° FOV. Maps due to the regular components 
in the model ADON (leftmost panel) and due to the random components in one of 200 
simulations for ADON30 (second panel) are shown. The combined map due to the regular 
and random components (third panel) as well as the binned map of median RMs in 2° x 2° 
pixels (rightmost panel) are also shown. We can see that the regular components of density 
and magnetic field alone produce large-scale structures of up to the FOV size. The structures 
highlight -By as induced by the radial and azimuthal components of the spiral and toroidal 
fields (the model shown does not include a poloidal field). The random components of density 
and magnetic field, on the other hand, produce complex structures, such as clump-like and 
filament-like features, on angular scales of ~ 1 — 10°. Such structures persist in the combined 
map. 



Mao et al.l (l2010l ) showed binned maps of median RMs in 2° x 2° pixels toward the NGP 
and SGP, which were produced from observational data. The maps display structures of 
a few to several degrees. In our model, such structures are mostly the consequence of the 
random components of density and magnetic field (rightmost panel). Our results indicate 
that the random components mainly produce the observed RM structures on scales less than 
several degrees, while the regular components contribute to larger scale structures. 



3.2. Contribution from Regular Components 

We first examine the statistics of RMs in maps due to the regular components alone. 
We calculated the average, fi = ^RM/Npi^, and the standard deviation, a = {^(RM — 
/i)^/(A^pix — 1)}^^^, in maps for the models and FOVs we consider. The resulting values are 
shown in Table 3. The average is mainly determined by the existence of the poloidal compo- 
nent of the GMF; /i ~ rad m~^ without the poloidal component, or /i ~ —4.8 rad m~^ and 
+5.6 rad m~^ toward the NGP and SGP, respectively, if the dipole poloidal component ex- 
ists, and fi ~ —4.1 rad m~^ and +4.7 rad m~^ toward the NGP and SGP, respectively, if the 
X-field poloidal component exists. The absolute values toward the NGP are a little smaller 
than those toward the SGP, due to the existence of a low density, local hot bubble centered 
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in the northern sky (Figure [2]). The averages for the 200 deg^ FOV are shghtly smaller than 
those for the 900 deg^ FOV, since the smaller FOV includes high latitude regions only where 
the density and magnetic field strength are both smaller. 

The standard deviations are a ~ 0.2 — 0.9 rad m~^ and ~ 0.7 — 4.3 rad m~^ toward the 
NGP and SGP, respectively, for the 900 deg^ FOV. These are caused mostly by contributions 
from the radial and azimuthal components of the GMF (see also the leftmost panel of Figure 
E]). The DS models tend to have larger standard deviations than the DN and QS models, 
because the DS models have larger Bq at mid and high altitudes as noted in Section 2.1. 
Larger values toward the SGP than those toward the NGP are again due to the existence of a 
local hot bubble near the Sun. The standard deviations for the 200 deg^ FOV are smaller by 
a factor of ~ 2 than those for the 900 deg^ FOV, since the smaller FOV includes a narrower 
range of RM values. 



3.3. Probability Distribution Function, Average, and Standard Deviation 



We next examine the statistics of RMs due to both the regular and random components. 
Figure [6] shows the probability distribution functions (PDFs) of RMs in maps of 900 deg^ 
FOV, averaged over 200 RM maps, for all m odels we conside r. The PDFs roughly follow 
the Gaussian dist ribution, as point ed out by IWu et al.l ( 120091 ). The figure also shows the 
observed PDFs of iMao et al.l (120101 ). The average and standard deviation of /i (the average 
of RMs for a map) and a (the standard deviation of RMs for a map) over 200 maps were 
calculated. Figure [7| shows the resulting values for 900 deg^ FOV as well as for 200 deg^ 
FOV. 

The value of /i, which approximates the peak positions for the nearly-symmetric PDFs 
in Figure O is determined mostly by the regular components. So the average of is close 
to yU of the regular components alone in Section 3.1; (/i) ~ rad m~^, without the poloidal 
component of the GMF, or (/i) ~ —5 rad m~^ and ~ +6 rad m~^ toward the NGP and 
SGP, respectively, if the poloidal component exists. As noted in Section 1, the medians of 
observed RMs toward the NGP and SGP are ~ 0.0 ±0.5 rad m~^ and ~ +6.3 ±0.5 rad m~^, 
respectively (IMao et al. I boioh . Hence, the models with poloidal components better reproduce 
the observed average toward the SGP, but the models without poloidal components are 
preferred for the NGP. On the other hand, the fluctuation of /i over the 200 maps is rather 
small; the standard deviations of /i for 200 maps are < 0.75 rad m~^ for models with Vrras = 
15 km s~^ and < 1.2 rad m~^ for models with V^ms = 30 km s~^, as shown with error 
bars in Figure [3 These values are much smaller than the difference between the observed 
averages, ~ 6.3 rad m~^, toward the NGP and SGP. This means that the difference cannot 
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be explained by the statistical fluctuation caused by the random components in our models. 
These results indicate that none of our models can simultaneously reproduce the observed 
averages of RMs toward the NGP and SGP. 

In Figure |6l the simulated PDFs have narrower and more sharply peaked proflles than 
the observed PDFs. The width of the PDFs is quantifled by a. In Figure [TJ the averages 
of (J are 0.9 — 1.3 rad m~^ toward the NGP and 1.3 — 4.3 rad m~^ toward the SGP for 
models with V^ms = 15 km s~^, and 1.7 — 2.1 rad m~^ toward the NGP and 2.4 — 4.8 rad m~^ 
toward the SGP for models with V^ms = 30 km s~^. The fluctuation of a in 200 maps is 
again small; the standard deviations of a for 200 maps are < 0.45 rad m~^ for models with 



15 km s ^ and < 0.81 rad m ^ for models with K 



30 km s ^, as shown with 



error bars. On the other hand, the estimations of iMao et al.l (|2010[ ) with their observed RMs 
are a ^ 9.2 rad m~^ and 8.8 rad m~^ toward the NGP and SGP. So our estimations of a 



with simulated RMs are substantially smaller than the values of iMao et al.l ( l2010l ). We argue 
that all of our models fail to reproduce not only the peak positions but also the widths of 
the PDFs of observed RMs toward the NGP and SGP (see further discussions in Sections 4 
and 5). 

From Figure [71 we see the average and standard deviation of /i are not sensitive to the 
size of FOV. But the values of a for the 200 deg^ FOV are somewhat smaller than those for 
the 900 deg^ FOV, as expected. 



3.4. Power Spectrum and Structure Function 

We also calculated the two-dimensional power spectrum (PS) and structure function 
(SF) with the sky map of RM. They tell us at which angular scales most power of RM 
resides and the spatial structure of RM decorrelates. 

Figure E] shows the PS of RMs due to both the regular and random components in maps 
of 900 deg^ FOV, averaged over 200 RM maps, for all models we consider. The PS toward 
the SGP are larger than those toward the NGP. This is again attributed to the presence of a 
local hot bubble in the northern sky. Toward the SGP, the PS for the models with a dipole 
toroidal fleld (D) are larger than those for the models with quadrupole toroidal field (Q). 
This is because the dipole models have larger Bq, as shown in Figure |31 (The dipole and 
quadrupole models are identical toward the NGP, as noted in Section 2.1.) The amplitude 
of the PS depends on the rms flow speed; the PS for V^ms = 30 km s~^ are larger by up to a 
factor of ~ 3 than those for l^rms = 15 km s~^. The slope also depends on the rms flow speed. 
Larger Kms's result in shallower proflles. This is because the density and magnetic fleld PS 
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of supersonic flows (VJ-ms = 30 km/s corresponds to Mrms ~ 1 — 2) are shallo wer than those 



of su bsonic flows (Kms = 15 km/s corresponds to Mrms ~ 0.5 — 1) (e.g., see iPadoan et al. 



20041 ). We quantified the slope of PS, a, defined with P{k) oc /c" in the angular range of 
~ 1 — 10°, and show the resulting values in Figure [HI for 900 deg^ FOV as well as for 200 deg^ 
FOV. These are the averages for 200 maps. The Figure confirms the dependence of a on 
Vrms] ct's for Vrras = 30 km/s are larger by up to ~ 0.4 than those for Vrms = 15 km/s for 
the models shown in Figure [8] (with 900 deg^ FOV). Some of the PS have slopes consistent 
with the Kolmogorov slope —5/3, but others have steeper slopes, indicating that the GMF 
model influences the slope of the PS. 

From the observer's point of view, the SF of RMs is a statistical quantity that is easier 
to obtain than the PS. The n-th order SF is defined as 

Sn{r) = {\RM{x+r) - RM{x)\')x (4) 

with r = I r| , where the subscript indicates the averaging over the data domain of x. Figure 
[To] shows the second-order SFs (n = 2) in maps of 900 deg^ FOV, averaged over 200 RM 
maps, for all models we consider. The SFs monotonically increase with angular separation, 
and reach up to ~ 10 rad^ m""^ toward the NGP and ~ several xlO rad^ m~^ toward the 
SGP, respectively, at the angular separation of ~ 10°. The SFs toward the SGP are larger 
than those toward the NGP. Toward the SGP, the SFs for the models with dipole toroidal 
fields are larger than those for the models with quadrupole toroidal fields. As for the PS, 
the SFs for V^ms = 30 km s~^ have amplitudes larger by up to a factor of ~ 3 than those for 
K-ms = 15 km s~"^. 



Figure do] also shows the observed SFs toward the NGP and SGP (IStil et al.ll201lf ). It 
is clear that the SFs of our simulated RMs do not match the observed SFs at any angular 
scales. At the largest observed scale of ~ 10°, the simulated SFs are smaller by an order of 
amplitude or so than the observed SFs. But the difference is even larger at smaller angular 
separations. Not just the amplitudes, but also the slopes of the SFs of our simulated RMs are 
quite different from the observed ones; the slopes of simulated RMs are much steeper. We 
quantify the slope of SF, defined with 5*2 (r) oc over the angular separations of ~ 1 — 10°, 
and show the resulting values in Figure for a 900 deg^ FOV as well as for a 200 deg^ FOV. 
Again these are the averages for 200 maps. For the 900 deg^ FOV, C, — 0.6 — 0.85 for 
Kms = 15 km s~^, and C ~ 0. 4 — 0.8 for V^ms = 30 km s^^. The slopes of observed SFs are 



~ 0.02 — 0.05 ( IStil et al.ll201l[ l. So our models fail to reproduce the observed SFs not only in 
the amplitude but also in the slope. The results suggest that the observed RMs may contain 
contributions due to structures smaller than those typically found in the Galactic halo. 

The PS and SFs for the 200 deg^ FOV (not shown) show behaviors similar to those for 
the 900 deg^ FOV, except that there are no data beyond a scale of ~ 7°. Quantitatively, 
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the slopes of PS and SF, a and C, for the 200 deg^ FOV are a bit smaller than those for the 
900 deg^ FOV, as shown in Figure |9l 



Discussion 



In our models, the average or median of the RM toward high Galactic latitudes, /i, is 
determined mostly by the halo poloidal component of the GMF, while the standard deviation, 
0", is determined mostly by the random component of the GMF. 

Our models fail to simultaneously reproduce the median s of observed RM s toward the 
NGP and SGP, ~ 0.0±0.5 rad m'^ and ~ +6.3±0.5 rad m^^ jMao et aliboiol ). The models 
that contain a poloidal component of the GMF have a vertical magnetic field of |-B|| | ~ 0.3 /xG 
in the Earth vicinity, which induces (/i) ~ —5 ± (0.3 — 0.5) rad m~^ toward the NGP and 
(yu) ^ 6 ± (0.6 — 1.1) rad m~^ toward the SGP. The models without a poloidal component 
have ':=L /xG, so (/i) ^ ± (0.4 — 1.2) rad m~^. The different values, |-B||| ^ /iG toward 
the NGP and ~ 0.3 /xG toward the SGP, which would explain the observed medians, however, 
is not accommodated in our r nodels. The diffic ulties of mixed regular field geometries in the 
steady state are discussed by iMao et al.l (|2010( ). 



The fiuctuation (standard deviation) of \x in simulations, < 1.2 rad m~^, is too small to 
account for the difference between the observed medians toward the NGP and SGP. To see 
whether larger fluctuations of /i are possible with different model parameters, simulations 
with larger values of Ldrive and Kms were carried out (not shown). We found that the 
fluctuation of /i increases with increasing Ljrive and becomes ~ 5 rad m~^ for Ldrive ^ 1-5 — 
2 kpc. However, with such large value of Ldrive, the EM distribution estimated from WHAM 
data (IHill et al.l 120081 ) cannot be explained (see Appendix). The fluctuation of /i increases 
with increasing V^ras also. We flnd a sufficiently large fiuctuation for Vrras ^ 100 km s~^, 
which corresponds to Mrms ~ 6 — 8. But such large Mrms is clearly inconsistent with the 
value constrained from the EM distribution (IHill et al.ll2008l ) and from polarization gradients 
JOaensler et al.ll201ll : ISurkhart et al1l2012h . 



Our models also fail to reproduce the standard deviations of obs erved RMs, a ^ 
9.2 rad m'^ and 8.8 rad m'^ toward the NGP and SGP, respectively jMao et allboioh . 
Our estimates, (a) < 4.5 ± (0.2 — 0.8) rad m~^, are substantially sr naller. The value s 
are not consistent even with the Galactic contributions estimated by ISchnitzelerl (120101 ). 
cr ~ 6.8 ± 0.1 rad m~^ and 8. 4 ifc 0.1 rad toward t he NG P and SGP, respectively. (Note 
that the data of iMao et al.l (l2010l ) and ISchnitzelerl (120101 ) are different.) We would get 
(cr) ~ 10 rad m~^ if we adopt V„ns ~ 100 km s~^. But again such a value is too large to 
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explain the observed EM distribution. 



We could explain the observed medians and standard deviations if the re are additional 
global or l ocal s tructures not represented in our models for the GMF (e.g., iMao et al.ll2010l : 



Stil et al.l I2OIII ). Such structures could have bee n produced, for instan ce, by the Parker 
instability or supernova explosions. In addition, iMachida et al.l ( 12012| ) have shown that 
transient structures can be produced in three-dimensional MHD simulations of turbulent 
gaseous disks. We could also address the discrepancy between the simulations and the 
data if the regular field is stronger than that in our models. For instance, we would have 
simulated the observed ex's, if the regular field strength was increased by a factor of ~ 
2. But then, the corresponding regular field strength near the Galactic plane would be 
4 (see Figure B]), which is larger than that adopted in re cent observational studies 



( Pshirkov et al. 2011 : Van Eck et al. 2011 : Jansson fc Farrai 2OI2I). O n the other hand, after 



investigating RMs toward the Perseus arm region, iMao et al.l (120121 ) argued that observed 
RMs are consistent with a toroidal field of strength ~ 2 yuG toward the north, but ~ 7 fiG 
toward the south (note that the observation covers a rather thin region of 8.8 kpc < R < 
10.3 kpc and 0.8 kpc < \z\ < 2.0 kpc). This indicates that there may be room for additional 
modeling of the GMF, which we leave for future studies. 

Even with additional s tructures and st ronger regular fields, however, the observed SFs 
are difficult to reproduce (IStil et al.l 1201 ll ). The second-order SFs of simulated RMs are 
substantially smaller than the observed SFs, especially at small angular separations of 0°.l — 
1°. This means that we need structures that would provide significant power s on scales of a 



parsec or so, but there is no observational support for such structures f e.g.. iGaensler et al. 



2005 



Haverkorn et al.ll2008D . The slope of observed SFs, 0.02 - 0.05 (IStil et al.ll201lD . is 
much fiatter than that of simulated SFs, ~ 0.4 — 0.9. We note that two-dimensional, white 
noise results in a fiat SF; 5*2 ~ 2(Tgj.j., where (Terr is the standard deviation of the noise. We 



200 rad^ m ^ for Up- ~ 



10 rad m ^. Observational errors should follow the 



would get 5*2 

distribution of white noise. But iMao et al.l ( 120101 ) estimated that the observati onal errors in 
their data are ~ 5 and ~ 3 rad m~^ toward the NGP and SGP, respectively. iTaylor et al. 



( 120091 ) estimated that the error in their data is ~ 8 rad m ^. So observation al errors may not 
be en ough to explain the observed fiat SFs of ~ 100 — 300 rad^ m~^ (see also lHammond et al. 
2012h . 



It is expected that radio sources against which RMs are observed have their own intrinsic 
RMs. The RMs of radio sources should also follow the distribution of white noise and 



contribute to a fiat SF. ISimonetti et al.l (119841 ) and ISimonetti fc Corded (119861 ) showed that 
the observed SF toward the NGP deviates from a fiat behavior at very small scales; the 
square value of RM difference between two sources at and Xj, {RM{xi) — RM{xj)}'^ (not 
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the second-order SF), is around ~ 10 rad^ at angular separations of ~ 0.01°. This tells 
us that even if the square value is solely due to the RMs of radio sourc es, it is unlikely that 
the RM s make up the observed flat SFs of ~ 100 — 300 rad^ m~^. In fact. lSimonetti fc Cordes 
(119861 ) claimed that many observed sources possess little intrinsic RM from the fact that the 
fractional polarization percentage is nearly constant with increasing wavelength. 

Finally, there exists clear evidences for an e xtragalactic component to the RMs, from 



three different pa pers taking di f ferent appr oaches (lSchnitzeleiil2010l : iBernet et al.ll2012l : iHammond et al, 
2012). Although Bernet et al. (2012) and Hammond et al. ( 20121 ) attribute the excess scat- 
ter in RM to individual intervening absorbers along the line of sight, an alternative origin 
which can contribute to the observed RM toward high Galactic latitudes too is the IGMF. 



As mentioned in Section 
universe 



, usin g simulations for the large-scale structure formation in the 
Akahori fc Rvul (l201lh predicted that the RM due to the IGMF would have /i ~ 



and a ~ several rad m~^. They also predicted that the SF of the RM has a flat profile of 
100 — 200 rad^ m"*^ amplitude at angular separations larger than 0.2°; the SF is expected to 
decrease to the order of ~ 10 rad^ m~^ at ~ 0-01° fro m ~ 100 rad^ m~'^ at ~ 0.1°, consis- 



tent with the observation of ISimonetti fc Corded (119861 ). The work of I Akahori fc Ryul (120111 ) 
suggests the possibility that a substantial fraction of the RM toward the NGP and SGP 
can be attributed to the RM due to the IGMF. If so, the observed standard deviations and 
SFs can be explained. But reproducing the observed medians of RM still needs additional 
components or/and structures of the GMF. 



5. Conclusion 

We have studied the Galactic RM toward high Galactic latitudes. We have considered 
a number of models for the global, regular components of the GMF and the electron density 
in the Galaxy, based on observations. The turbulent, random components were modeled 
with three-dimensional MHD turbulence simulations. The strength of the regular magnetic 
field in our models is a few /xG close to the disk and smaller at high altitudes. The strength 
of the turbulent field is at most a few /iG, and mostly < 1 /iG. We obtained RM maps for 
900 deg^ FOV toward the Galactic poles, and compared the results with observations. We 
also considered a smaller FOV of 200 deg^, designed to simulate the FOV of future surveys 
with the SKA. 

Our models fail to simultaneously reproduce the observed medians of RMs toward the 
NGP and SGP. The observations require vertical magnetic fields of B\\ ~ /uG toward the 
NGP and 0.3 /xG toward the SGP in the Earth's vicinity, but such field geometries are not 
accommodated in the GMF models we considered. The PDFs of simulated RMs are narrower 
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and more sharply peaked than the observed PDFs, meaning that the standard deviations 
of simulated RMs are smaller than the observed values. The second-order SFs of simulated 
RMs are one to two orders of magnitude smaller than the observed SFs at small angular 
separations. In addition, the slopes of the SFs of simulated RMs are substantially larger 
than the observed ones. 

We argue that observational errors and the intrinsic RM of background radio sources 
are not enough to explain the discrepancies between the statistics of our simulated RMs and 
observed RMs. We suggest that the RM due to the IGMF may account for a fraction of 
the RM toward high Galactic latitudes and could explain the discrepancies in the standard 
deviation and SF of the RMs. As a subsequent project, the RM due to the IGMF as well 
as observational errors and the intrinsic RM of background radio sources would be incor- 
porated in the modeling of the RM toward high Galactic latitudes, to check quantitatively 
whether some of the discrepancies, such as those in the standard deviation and SF, would be 
explained. But we expect that reproducing the observed medians will still need additional 
components or/and structures of the GMF. 

We should note, however, that current observations of RMs still contain large uncertain- 
ties. New and the future observational facilities such as the JVLA, MWA, LOFAR, ASKAP, 
MeerKAT, and the SKA will produce much better data. For instance, ASKAP will detect 
radio sources with average angular separations of ~ 0.1°. With such a dense RM grid, the 
quality of the observed RM data will be dramatically improved. Better quality data will 
hopefully enable us to better quantify the contributions due to the GMF, the IGMF, the 
intrinsic RM of background radio sources, and observational errors. 

Finally, we note that better measur ements of the Mach number of turbulence in the 



halo by observations of Hq, line profiles (IHaffner et al. 



2OIOI ) as well as by radio polarization gradients (IGaensler et al. 



2003 



Hi 



1 et al. 2008; Haffner et al. 



20111 ) will help us improve 



the constraints on the magnitude and structures of the RM toward high Galactic latitudes. 



A. Coordinate Systems 

The following coordinate systems are used throughout the paper: Cartesian coordinates, 
{x,y,z), Galactocentric cylindrical coordinates, {R,Q,z), Galactocentric polar coordinates, 
{r,9,(l)), and Galactic celestial coordinates, {l,b), respectively, defined in Figured! Here, 
R = (x^ + y'^Y^'^ is the Galactocentric radius, = is the azimuth angle starting from 
/ = 90° and increasing in the counterclockwise direction, and the x-y plane coincides with 
the Galactic plane with x pointing to I = 90° and y to I = 180°. The Sun is located at 
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(x, y, z) = (0, 8.5, 0) in kpc, and thus Rq = 8.5 kpc. 



B. Electron Density 



The NE2001 model we employ for the electron density, ne,o, consists of three global 
components (thick disk, thin disk, and five spiral arms), four local components (local hot 
bubble. Loop 1, local super bubble, and low density region in quadrant 1), and more than 
a hundred isolated components such as clumps and voids, as displayed in Figure [TJ The 
electron density of the model was obtained from the average of the DM for pulsars at known 
distances. We take the original parameters in the NE2001 packag^, except for the scale 
height of the thick d isk, hi, and the mid-plane electron density for the thick disk, ne,i- 
Gaensler et al.l (120081 ) obtained hi ~ 1.8 kpc from an analysis of pulsars at high Galactic 
latitudes. We adopt hi = 1.8 kpc, instead of 0.97 kpc in the package. The corresponding 
mid-plane electron density for the thick disk is ne,i = 0.014 cm~^. 

The revised scale height improv es fits of the Galacti c RM and radio continuum emission 
at low and mid-Galactic latitudes (ISun &: ReichI l2010l ). The changes of the scale height 
and mid-plane electron density of the thick disk, however, may not be consistent with the 
construction of the NE2001 model, because these parameters are highly covariant with others 
in the model. In order to justify our modification and make sure the consistency with 
observations, we compare the DM and EM from the NE2001 model with the original scale 
height and mid-plane density (hereafter the original NE2001 model) and the NE2001 model 
with the modified scale height and mid- plane density (our mo dified NE2001 model). We also 
tested the exponential model fitted by iGaensler et al.l (120081 ) (the plane-parallel model) for 
comparison. 

Figure [TT] shows the distribution of DM as a function of height above the Galactic plane. 
At a height of ( 3.2 — 0.4 kpc, both th e original and modified NE2001 models reproduce the 
distribution of iGaensler et al.l ( 120081 ) well toward the SGP, but predict smaller DM by a 
factor of ~ 2 toward the NGP, likely due to the local hot bubble. At the height of > 1 kpc, 
the original NE2001 model overestimates DM by a factor of ~ 1.5 toward south. On the other 
hand, the modified NE2001 model reproduces the observed DM well; the resultant DMs for 
the modified NE2001 are 22 pc cm~^ and 25 pc c m~^ toward the NGP and SGP, respectively, 



which are roughly c onsistent with observations (jPeterson fc Webbeii l2002l : iHill et al.l 12008 
Gaensler et al]l2008h . 



^We use NE2001_1.0 downloaded from 



http://astrosun2.astro.cornell.edu/~cordes/NE2001/ 
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Next, we calculated the PDF of EM. There is a well-known issue that models with 
smooth electron density profile do not reproduce the observed EM well; the discrepancy can 



be resolved by introducing volume filling factor, /e. We adopted fp = 0.07 exp 
in kp c, for z < 0.75 kpc and /e = 0.32 for z > 0.75 kpc (iBerkhuijsen et al.l 12006 



Sun et al. 



20081). The resulting EM is shown in Figure [H With the adopted /e, the origin al NE2001 



does not correctly reproduce the EM distribution from the WHAM observation ( iHill et al. 



20081 ). On the other hand, the modified NE2001 model and the plane-parallel model produce 



results, which are better consistent with the observed distribution. 

It is interesting to see that the cases with V„ns = 15 km s~^ (or the rms Mach number 
of t urbulence is clo s e to ^ 1) better reproduce the observed EM. It should be pointed 
that iGaensler et al.l (|2008[ ) claimed smaller /g (that is, more clumpy) at both low and high 
altitudes. If smaller /g is adopted, our model gives larger dispersion of EM, so that even 
smaller V^rms would be preferable. Then, the corresponding standard deviation of RM would 
be even smaller (Figure [7]). 

Finally, Figure [12] shows the second-order SF of simulated RMs in a 900 deg^ FOV. 
The amplitude of SF for the original NE2001 model is somewhat larger than that for the 
modified NE2001 and plane-parallel models. The slope of SF does not significantly depend 
on the electron density model, and again the SF for the three models does not match the 
observed one (see Section 3.4). 



C. Magnetic Field 



The global, regular magnetic field, Bq, is conventionally described as the combination 
of the disk spiral field, Bg, the halo toroidal field, and the halo poloidal field. Bp, so that 



Bs + Bt + Bp (e.g., iProuza fc Smidall2003l ). There are a number of analytic models 



Bo 

intended to reproduce the observed structure of Bq. The models have been tuned by fitting 
the data mostly at low and mid-Galactic latitudes, but can be used for studies of the Galactic 
R M toward h i gh Ga lactic latitudes as well. Specifically, we emp loy the models introduced 
bv ISun et al.l ( 120081) for the spiral field, by ISun fc ReichI (l2010r) for the toroidal field, by 



Giacinti et al.l (120101 ) for the dipole poloidal field, and by IJansson fc Farrarl (120 12[ ) for the 
X-field poloidal field, unless otherwise specified. We note that our combined model does 
not take account of detailed theoretical consistencies, such as the continuity, divergence-free 
nature, and closeness of Bq. We expect that those do not significantly affect our results. 

The large-scale disk fields of spiral galaxies including our Galaxy have been classified 
into two types based on the spiral pattern, the axi-symmetric spiral (ASS) field with no 



-20- 



dependence on the azimuthal angle, or the bi-symmetric spiral (BSS) field with a symmetry 
of 7T. The functional form for the disk spiral field can be written as 



B,^r{R, e, z) = Di{R, z)D2{R, 9) sin(pso), 
5,,e(i?,e,^) = -D,{R,z)D2iR,Q)cosip,o) 
B,,,{R,Q,z) =0, 



(CI) 



where Pso is the pitch angle of arms which is positive for leading spirals and negative for 
trailing spirals such as in our Galaxy, and 



Di{R,z) 



Bso exp 
Bsc exp 



— 5 — - — — ) R> Rsc, 



(C2) 



\z\ 



R < Rfic 



For the ASS model, we adopt the ASS + reversals of the magnetic field directions: 



D2{R,Q) = I 



+ 1 R> 7.5 kpc, 

-1 6 kpc < i? < 7.5 kpc, 

+1 5 kpc < i? < 6 kpc, 

-1 R<5 kpc. 



(C3) 



where +1 means the clockwise direction as seen from the north pole. We also adopt Rgo 
10 kpc, Zso = 1 kpc, -Rsc = 5 kpc, -Bso = 2 /iG, Bsc = 2 fiG, and pso = —12°. 



For the BSS model, we adopt the spiral structure of 

D^iR, 6) = sin (q + In ^ 

V tanpso ^sb 



(C4) 



We also adopt Rso = 6 kpc, Zso = 1 kpc, Rsc = 3 kpc, Bso = 2 /iG, and Bsc = 2 fiG. 
In addition, Rsh = 9 kpc and pso = —10° for i? > 6 kpc, and otherwise Rsh = 6 kpc and 
Pso = —15°. Note that the trigonomet ric function is ra inus cosine if the azimuth angle, G, 
starting from I = 180° is adopted (e.g.. lSun et al.ll2008[ ). 



Reversals in th e sign of RM across the Galactic plane and across the Galactic center (e.g., 
Taylor et al.l 120091 ) suggest the existence of a halo t oroidal field. Global three-dimensional 
MHD simulations of gas disks (iNishikori et al.ll2006[ l have indicated the azimuthal magnetic 
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field component in the halo as a result of the buoyant escape of the azimuthal magnetic flux 
from the disk. The halo toroidal (azimuthal) field can be expressed as 



BtAR, e, z) 

5t,e(i?,e,z) 



0, 

signjzyBto R 

I _^ f \z\-zto \'^ Rto 

0, 



exp 



R-R 



to 



R 



to 



(C5) 



where sign(z) is the sign of z, and v introduces the parity of the toroidal field configuration; 
w = 1 if we consider the asymmetries in longitude and latitude relative to the Galactic plane 
and the center, respectively (dipole), or f = 2 if we consider the axisymmetric configuration 
without reversals relative to the Galactic plane (quadrupole). We adopt Bto = 2 /xG, Zto = 
1.5 kpc, Rtn = 4 kpc, and_Zti = 0.2 kpc for \z\ < Zto and zu = 4.0 kpc otherwise. We note 



that IProuza fc Smidal ( 20031 ) adopted smaller zn for high z (and also similarly small Zt2 in 
their mo del), which results in Bt of order 0.1 fiG at z ~ 3 kpc. We adopt the parameters 
based on ISun fc ReichI (120101 ) . which gives ~ 1 fiG at high z (see also Figure [3]). Such a 
strong halo field at high z is also motiv ated by recent studies of halo magnetic fields (see 
Mao et all [201a : Ijansson fc Farrai]l2012h . 



Nonthermal filamen ts observed near the Galactic center imply the existence of a central, 
vertical magnetic field (jHan 2009 ). which is also predic ted from global three-dimensional 



MHD simulations of gas disks (e.g., iMachida et al. 



2OO9I ). The vertical field can be due to 



the Galactic-center poloidal (dipole) field; a strong dipole field would be observed as the 
vertical magnetic field in the Earth vicinity. The dipole field can be expressed as 



B^^R{R,e,z) 



3Rz 



(i?2 + z2)3/2 R2 + 



' Bp,e{R,Q,z) =0, 



5p,.(i?,e,^) 



(C6) 



(/?2 + ^2)3/2 yR2^^2 



3z^ 



Here, we set fip = 180 fiG ■ kpc^ (iGiacinti et al.l 120101 ) to make the vertical component of 
~ 0.3 fiG from the poloidal field. The verti cal field could be due to the out-of-plane X-field 
recently studied by lJansson fc Farraii (120121 ). This model is motiv ated by the X-shaped fiel d 
structure seen in radio observations of external, edge-on galaxies (lBeckll2009l : lKrausell2009l ). 
We define the elevation angle of magnetic field, r]{R,Q,z), with respect to the Galactic 
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mid-plane. Then, the X-field can be expressed as 




w 



B^,R{R,e,z) 



sign{z)B^oexp 



cos 77, 



(C7) 



( 




w 



B^,,{R,e,z) 



-Bxo exp 



sin?7. 



We take the field at Rp > R^c to have a constant elevation angle, rj^o, with respect to the 
mid-plane, where _Rp is the radius at which the field line passing through {R, 0, z) crosses 
the mid-plane. This means that if R—\z\/ tan?7xo > -RxO) then Rp = R—\z\/ tan?7xo, f] = rj^o, 
and w = 1. Otherwise, Rp = RR^c/iRxc + \AI tanr^xo) 'f] = tan^^{|z|/(i? — Rp)} and w — 2. 
We adopt B^o = 4.6 /xG, rj^o = 49°, i?xc = 4.8 kpc, and i?xo = 2.9 kpc. 
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Table 1. Models, regular field direction, Mach number, and 



ModeP Block 1 Block 2 Block 3 Block 4 

(0.00 < \z\ < l.2by= (1.25 < \z\ < 2.50)^ (2.50 < l^l < 5.00)^ (5.00 < \z\ < 10.0)"° 



AD0N15 


+.T, 


1.0, 


3.0 


— X 


1.0, 


3.0 


— X 


0.5, 





1 


— X 


0.5, 


0.1 


AD0S15 


+X, 


1.0, 


1.0 


+x, 


1.0, 


0.1 




0.5, 





1 


+x, 


0.5, 


0.1 


AQ0S15 


Jf-X, 


1.0, 


3.0 


— X 


1.0, 


3.0 


— X, 


0.5, 





1 


— X, 


0.5, 


0.1 


ADPN15 


+r. 


1.0, 


1.0 


— X 


1.0, 


1.0 


— X 


0.5, 





1 


— X 


0.5, 


0.1 


ADPS15 


+.T. 


1.0, 


1.0 


+x, 


1.0, 


0.1 


+x, 


0.5, 





1 


+x, 


0.5, 


0.1 


AOPS15 


+ T 


1.0, 


1.0 


— X 


1.0, 


1.0 


— X 


0.5, 





1 


— X 


0.5, 


0.1 


ADXN15 


+x, 


1.0, 


1.0 


— X. 


1.0, 


1.0 


— X, 


0.5, 





1 


— X, 


0.5, 


0.1 


ADXS15 


+x, 


1.0, 


0.1 




1.0, 


0.1 


+x, 


0.5, 





1 


+a;, 


0.5, 


0.1 


AQXS15 


+x. 


1.0, 


1.0 


— X 


1.0, 


1.0 


— X 


0.5, 





1 


— X. 


0.5, 


0.1 


BD0N15 


— X, 


1.0, 


3.0 


— X, 


1.0, 


3.0 


— X, 


0.5, 





1 


— X, 


0.5, 


0.1 


BD0S15 


+x, 


1.0, 


1.0 


-\-x, 


1.0, 


0.1 




0.5, 





1 


+a;, 


0.5, 


0.1 


BQ0S15 


— X 


1.0, 


3.0 


— X 


1.0, 


3.0 


— X. 


0.5, 





1 


— X 


0.5, 


0.1 


BDPN15 


— X, 


1.0, 


1.0 


— X 


1.0, 


1.0 


— X. 


0.5, 





1 


— X. 


0.5, 


0.1 


BDPS15 


+x. 


1.0, 


1.0 


+x, 


1.0, 


0.1 


+x, 


0.5, 





1 


-\~x, 


0.5, 


0.1 


B0PS15 


— X 


1.0, 


1.0 


— X 


1.0, 


1.0 


— X 


0.5, 





1 


— X 


0.5, 


0.1 


BDXN15 


— X. 


1.0. 


1.0 


— X 


1.0, 


1.0 


— X 


0.5, 





1 


— X 


0.5, 


0.1 


BDXS15 


-\-x, 


1.0, 


1.0 


-\-x, 


1.0, 


0.1 


-\-x, 


0.5, 





1 


+x, 


0.5, 


0.1 


BQXS15 


— Xf 


1.0, 


1.0 


— X. 


1.0, 


1.0 


— X, 


0.5, 





1 


— X, 


0.5, 


0.1 


ADON30 


+x, 


2.0, 


1.0 


— X, 


2.0, 


1.0 


— X 


1.0, 





1 


— X, 


0.5, 


0.1 


ADOS30 


+x, 


2.0, 


1.0 


+x, 


2.0, 


0.1 


+x, 


1.0, 





1 


+x, 


0.5, 


0.1 


AQOS30 


+x, 


2.0, 


1.0 


—x, 


2.0, 


1.0 


— X, 


1.0, 





1 


— X, 


0.5, 


0.1 


ADPN30 


+x, 


2.0, 


1.0 


—X, 


2.0, 


1.0 


— X, 


1.0, 





1 


— X, 


0.5, 


0.1 


ADPS30 


+x, 


2.0, 


1.0 


+x, 


2.0, 


0.1 


+x, 


1.0. 





1 


+x, 


0.5, 


0.1 


AQPS30 


+x, 


2.0, 


1.0 


—X, 


2.0, 


1.0 


— X, 


1.0, 





1 


— X, 
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ADXN30 
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— X, 


1.0, 





1 


— X, 
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0.1 


ADXS30 
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2.0, 


0.1 


+X, 


2.0, 


0.1 


+x, 


1.0, 





1 
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0.5, 


0.1 


AQXS30 


+x, 
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1.0 


—X, 


2.0, 


1.0 


— X, 


1.0. 





1 


— X, 


0.5, 


0.1 


BDON30 


— X, 


2.0, 


1.0 


—X, 


2.0, 


1.0 


— X, 


1.0, 





1 


— X, 


0.5, 


0.1 


BDOS30 


+x, 


2.0, 


1.0 


+x, 


2.0, 


0.1 


+x, 


1.0, 





1 


+x, 


0.5, 


0.1 


BQOS30 


— X, 


2.0, 


1.0 
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1.0 
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1.0, 
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0.5, 


0.1 


BDPN30 


— .T, 
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— X, 
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1 
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0.5, 


0.1 
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~x, 
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—X, 


2.0, 
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1 
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0.1 
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1.0 


+x, 
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0.1 


+x, 
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1 


+x, 
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Table 1 — Continued 



ModeP 


Block 1 


Block 2 


Block 3 


Block 4 




(0.00 < l^l < 1.25)'' 


(1.25 < 1^1 < 2.50)'' 


(2.50 < \z\ < 5.00)^ 


(5.00 < \z\ < 10.0)^ 




-Bo, Mrma, 00 




Bo, Mj-jns, A 


Bo, Mrms, Po 


BQXS30 


-X, 2.0, 1.0 


-X, 2.0, 1.0 


-X, 1.0, 0.1 


-X, 0.5, 0.1 



''^A: axi-symmctric spiral, B: bi-symmetric spiral, D: dipole toroidal, Q: quadrupole toroidal, O: 
no poloidal, P: dipole poloidal, X: X-field poloidal, 15: random field with rms flow speed 15 km/s, 
30: random field with rms flow speed 30 km/s, N: toward the NGP, S: toward the SGP. 



''In units of kpc. 
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Table 2. Rms values of density and magnetic field strength in MHD turbulence 

simulations. 





/3o 


Prms/ PO 


-Brms/ -Bo 


0.5 


0.1 


1.02702 


1.00360 




1.0 


1.03139 


1.03748 




10. 


1.02755 


1.34981 


1.0 


0.1 


1.09307 


1.01275 




1.0 


1.11417 


1.14356 




3.0 


1.10303 


1.40630 




10. 


1.11177 


2.09291 


2.0 


0.1 


1.29744 


1.04641 




1.0 


1.32886 


1.45009 




10. 


1.32293 


3.04786 
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Table 3. The average, and standard deviation, cr, of the RMs due to the regular 
components of the electron density and the GMF. 





200 deg^ 


900 deg^ 


Mndpl 


jj,^ 




IJ,^ 


a'' 


ADON 


0.00 


0.18 


0.06 


0.44 


ADOS 


0.02 


1.92 


0.03 


4.21 


AQOS 


0.02 


0.39 


0.03 


0.90 


ADPN 


-4.76 


0.45 


-4.83 


0.87 


ADPS 


5.58 


1.95 


5.66 


4.22 


AQPS 


5.58 


0.36 


5.65 


0.78 


ADXN 


-4.12 


0.13 


-4.12 


0.29 


ADXS 


4.69 


1.97 


4.70 


4.28 


AQXS 


4.69 


0.51 


4.70 


1.11 


BDON 


0.00 


0.14 


0.05 


0.35 


BDOS 


0.02 


1.87 


0.03 


4.08 


BQOS 


0.02 


0.32 


0.03 


0.75 


BDPN 


-4.76 


0.41 


-4.83 


0.80 


BDPS 


5.58 


1.90 


5.65 


4.09 


BQPS 


5.58 


0.31 


5.65 


0.66 


BDXN 


-4.12 


0.10 


-4.12 


0.24 


BDXS 


4.69 


1.91 


4.70 


4.14 


BQXS 


4.68 


0.44 


4.70 


0.96 



^In units of rad m 




Fig. 1. — Electron density distribution, and definitions of tlie coordinate systems used in 
tfiis paper. Small boxes in the bottom-right panel depicts the configuration of the data 
stacking (Section 2.2.2). 
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Fig. 2. — One-dimensional profiles from the Sun along the Galactic poles. Shown are the 
regular electron density (top left, thin red and thick blue are toward the NGP and SGP, 
respectively), the electron temperature (bottom left), the rms speed of random flow motions 
(top right), and the rms Mach number (bottom right, thin and thick lines are for V„ns = 15 
and 30 km s~^, respectively). 
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Fig. 3. — One-dimensional profiles from the Sun toward the Galactic poles. The top panels 
show the regular field strength, Bq (solid), and the LOS field strength, B\\ (dashed). B\\ is 
pointing away from us and toward us for models toward the NGP and SGP, respectively. 
The bottom panels show the plasma beta of the regular magnetic field, /3o- Profiles only 
for axi-symmetric spiral models are shown. The names of models mean: (A) axi-symmetric 
spiral, (B) bi-symmetric spiral, (D) dipole toroidal, (Q) quadrupole toroidal, (O) no poloidal, 
(P) dipole poloidal, (X) X-field poloidal, (N) toward the NGP, and (S) toward the SGP. 




Fig. 4. — One-dimensional profiles from the Sun toward the north Galactic pole. The left 
and right panels show the profiles for ADONIS and ADON30, respectively. The Panels from 
top to bottom show the electron density, the total magnetic field strength, the LOS magnetic 
field strength, and the cumulative RM from the Sun, respectively. The arrows with red and 
blue colors indicate the direction of magnetic field and the sign of RM. In the top two panels, 
the black lines show the regular components. B\\ of the regular magnetic field is zero along 
the Galactic poles in ADONIS and ADON30 without the poloidal component. 
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Fig. 5.— Two-dimensional RM map of the NGP with 30° x 30° FOV for ADON30. Panels 
from left to right show RM maps due to: the regular components of the electron density and 
the GMF only; the random components only; both the regular and random components; 
and a map binned in 2° by 2° pixels. The Galactic celestial coordinates are shown in the 
leftmost panel. 
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Fig. 6. — Probability distribution functions (PDFs) of simulated RMs in 30° x 30° FOV for 
models we consider (thin lines). Shown PDFs are the averages for 200 maps, where error 
bars indicate the standard deviation. Also shown as thic k lines are the P DFs of observed 
RMs toward the NGP (red) and SGP (blue), respectively jMao et al.ll201o[ ). 
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Fig. 7. — Average fi and standard deviation a of simulated RMs in 30° x 30° (left) and 
14.14° X 14.14° (right) FOVs for models we consider. Symbols are the averages of /i and a 
for 200 maps, while error bars indicate the standard deviations over 200 maps. Also shown 
as gray solid and dashed lines are observed values and errors: fi ~ 0.0 ± 0.5 rad m~^ and 
cr ^ 9.2 r ad m~^ toward th e NGP and fi ~ +6.3 ± 0.5 rad m~^ and a ^ 8.8 rad m~^ toward 
the SGP jMao et al.lboioi ). 
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Fig. 8. — Power spectra (PS) of simulated RMs in 30° x 30° FOV for models we consider. 
PS shown are the average over 200 maps, where error bars indicate the standard deviation. 
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Fig. 9.— Slopes of PS (a) and SF (C) of simulated RMs over 30° x 30° (left) and 14.14° x 
14.14° (right) FOVs for models we consider. Symbols are the averages over 200 maps, while 
error bars indicate the corresponding standard deviations. 
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Fig. 10. — Second-order structure functions (SFs) of simulated RMs over a 30° x30° FOV for 
models we consider. Shown SFs are the averages for 200 maps, where error bars i ndicate the 



standa rd deviatio n for 200 rnaps. A lso shown are the observed second-order SFs (IMao et al. 



(l2010l . circles) and lStil et al.l ( 1201 ll . lines)). Open circles and thick lines are toward the NGP, 
and filled circles and thin lines are toward the SGP. 
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Fig. 11. — Distributions of DM as a function of height above the Galactic plane. Distri- 
butions shown are the averages over 200 maps, where error bars indicate the corresponding 
standard deviations. Left to right panels show the results for the original NE2001 model, the 
modified NE2001 model (our model), and the plane-parallel model, respective ly. Symbols de- 
note t he pulsar observations at high galactic latitude (40 - 90 degrees) used by lGaensler et al. 
(120081 ): filled and open circles indicate observations with distance determination from trigono 
metric parallaxes and from associations with globular clusters, respectively. 
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Fig. 12. — Probability distribution functions (PDFs) of EM over a 60° x 60° FOV toward 
the north and south Galactic poles. Shown PDFs are the averages for 200 maps, where error 
bars indicate the standard deviation. Left to right panels show the results for the original 
NE2001 model, the modified NE2001 model (ou r model), and the plane-p arallel model. We 
adopted the volume filling factor introduced by iBerkhuijsen et al.l ( 2006i). Filled ci rcles are 
the WHAM observations at high (60 - 90 degree) galactic latitude (IHill et al.lboOSi ). 



-43 - 



100 _ 



■>a- 

E 

fN 

"§ 10 




E2001 
h=0.97 kpc 

ADONIS - 
ADON30 
ADPS15 
ADPS30 



J I I I L_ 




ADONIS - 
ADON30 

— ADPS15 

— ADPS30 

J I I I L 




1 10 
r [degree] 



1 10 
r [degree] 



1 10 
r [degree] 



Fig. 13. — Second-order structure functions (SFs) of simulated RMs over a 30° x 30° FOV. 
Shown SFs are the averages for 200 maps, where error bars indicate the standard deviation. 
Left to right panels show the results for the original NE2001 model, the modified NE2001 
model (our model), and the plane-parallel model. Open circles and black thick lines are the 
observed second- order SFs toward the NGP, an d filled cir c les an d black thin lines are those 
toward the SGP jMao et all J2OI0I . circles) and lstil et all (|201l[ lines)). 



